The impact of dissipation and noise on fluctuations in chiral fluid dynamics 



(N 



in 
o 



>< 



Marlene Nahrgang, 1,2 Christoph Herold, 2,3 Stefan Leupold, 4 Igor Mishustin, 2 and Marcus Bleicher 2 ' 3 

1 SUBATECH, UMR 6457, Universite de Nantes, Ecole des Mines de Nantes, 
IN2P3/CRNS. 4 rue Alfred Kastler, 44307 Nantes cedex 3, France 
2 Frankfurt Institute for Advanced Studies (FIAS), 
Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany 

3 Institut fur Theoretische Physik, Goethe- Universitat, 
Max-von-Laue-Str. 1, 60438 Frankfurt am Mam, Germany 
^Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden 

(Dated: January 20, 2013) 

We investigate the nonequilibrium evolution of the sigma field coupled to a fluid dynamic expan- 
sion of a hot fireball to model the chiral phase transition in heavy-ion collisions. The dissipative 
processes and fluctuations are allowed under the assumption that the total energy of the coupled 
system is conserved. We use the linear sigma model with constituent quarks to investigate the 
effects of the chiral phase transition on the equilibration and excitation of the sigma modes. The 
quark fluid acts as a heat bath in local thermal equilibrium and the sigma field evolves according to 
a semiclassical stochastic Langevin equation of motion. The effects of supercooling and reheating of 
the fluid in a first order phase transition are observed via the delayed relaxation of the sigma field 
to a new equilibrium state. Nonequilibrium fluctuations of the sigma field in a scenario with a first 
order phase transition are much stronger than in a critical point scenario. 
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I. INTRODUCTION 

73 ' 

At high temperatures and densities strongly interacting matter is supposed to form a state called the quark-gluon 
i " i , plasma, which can experimentally be created and studied in relativistic heavy-ion collisions. A broad variety of 
experimental signatures is proposed to study the QCD phase diagram in the laboratory. These signatures cover 
different aspects of the deconfmement and the chiral phase transition, including a suggested critical point [l|, 0] and 
the first order phase transition [1,01 • It is, however, not clear how much of a potential signal of the phase transition 
is developed in a realistic scenario of a heavy-ion collision. The fast dynamics during the creation and expansion of 
the fireball of strongly interacting matter makes it difficult to defer results from fluid dynamic calculations directly. 
A thorough theoretical understanding of phase transitions in the environment of a heavy-ion collision is necessary to 
make profound predictions for upcoming experiments at the low energy RHIC program [5(, the CBM project at FAIR 
Q and at NICA Q experiments. At a critical point the relaxation time for an order parameter becomes infinite. 
Therefore, a system that cools rapidly through a second order phase transition is necessarily driven out of equilibrium 
and any signal for a critical point is weakened Q. A fast expansion, however, leads to interesting phenomena at 
a first order phase transition, such as supercooling Q and spinodal decomposition [Toj . An intersting signal was 
proposed for a nonequilibrium chiral phase transition, i.e. the enhancement of soft pions from the decay of disoriented 
chiral condensates (DCC) [lll - ll3| . The theoretical description of such nonequilibrium processes can be achieved by a 
Langevin equation. It has been derived explicitly within a nonequilibrium quantum field formalism for the (p 4 model 
[I3 - ll7| . gauge theories [H, [l9| and O(N) chiral models [20]. Phenomenological treatments of stochastic Langevin 
equations for the propagation of the chiral fields have been performed to study the formation of disoriented chiral 
condensates [l2| or to study hadron production in a spinodal decomposition scenario [2l|. These studies were done 
at either a fixed temperature or assuming a scaling expansion and cooling of the system. 

In the following we address the two questions: how much of the equilibrium behavior at a critical point survives in 
a dynamical environment and how much of nonequilibrium effects at a first order phase transition survive relaxational 
dynamics? 

Our calculations below are performed within a chiral fluid dynamic model, where the propagation of chiral fields 
is coupled to a fluid dynamical evolution of the medium. In [22] this coupled dynamics was desribed within the 
linear sigma model with constituent quarks. In this model, the fluid consists of quarks and antiquarks which are 
considered in local thermal equilibrium and propagated according to energy-momentum conservation of an ideal fluid. 
Alternatively, in [23l - |25j the quarks were described by a Vlasov equation in the collisionless approximation. In (2(| [27j 
initial fluctuations of the sigma field were introduced and propagated deterministically during the evolution of the 
system. Then large oscillations of the chiral fields were observed, which did not show relaxational behavior of the 
order parameter. In the present work we go beyond deterministic equations of motion for the chiral fields by including 
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dissipation and noise in the framework of a Langevin equation. As already mentioned above we will utilize the linear 
sigma model with constituent quarks for this purpose. 

This paper is organized as follows: in section |n] we present the coupled model of extended chiral fluid dynamics 
starting from the linear sigma model with constituent quarks and in section IIIII we explain details of its numerical 
implementation with special emphasis on the energy-momentum conservation during the evolution of the entire system. 
In section HVl we investigate the nonequilibrium evolution of the coupled dynamics leading to supercooling, reheating 
and an enhanced intensity of the nonequilibrium fluctuations of the sigma field. Section [V] is devoted to summary 
and outlook. 



II. EXTENDED CHIRAL FLUID DYNAMICS 



A. General remarks 

Our work is an extension of existing chiral dynamic models [HI, [13, H2, Hfl including relaxational dynamics of 
the order parameter of the chiral phase transition in order to study the effects of nonequilibrium fluctuations. Due 
to the interaction with the quark fluid the chiral fields are damped. This is taken into account by friction terms in 
the Langevin equation of motion of the chiral fields. Due to these terms, the energy of the chiral field fluctuations 
dissipates into the quark fluid. To fulfil the dissipation-fluctuation theorem [28[ the field fluctuations are included as 
a noise term on the right hand side of the Langevin equation, which reads 

d^a+ 5 -^+ S -^+r,d t <T = Z. (I) 
ocr da 

It contains a classical Mexican hat potential U, the quark contribution to the thermodynamic potential VLq q to one- 
loop level, the damping coefficient r\ and the stochastic noise field £ 1 . The dynamics of the quarks is reduced to the 
propagation of densities according to energy-momentum conservation, i.e. the equations of relativistic fluid dynamics 

g^u = gu > (2) 

where the source term S v accounts for the energy-momentum exchange between the fluid and the field. 

Our model is a combination and an extension of both approaches: the chiral fluid dynamics and the Langevin 
description of the chiral fields at the phase transition mentioned above. This model was derived in [29| within the 
formalism of the two-particle irreducible (2PI) effective action [30h38j. which gives the Langevin equation for the 
sigma field and the dynamics of the quarks consistently. The field and the fluid are coupled at different levels of the 
evolution of the entire system: 

• The one-loop effective potential V e s = U + ilq q for the sigma field in presence of the quarks drives the chiral 
phase transition. Variation with respect to the sigma field gives the one-loop scalar density d g^ 9 = gp s , which 
is the standard mean-field contribution to the equation of motion ([T]) of the chiral field. 

• The pressure and the energy density of the locally equilibrated quark fluid in the equation of state p(e) depend 
on the local value of the sigma field, which plays the role of an external parameter in the thermodynamic sense. 

• The damping coefficient 77, which describes the dissipative and stochastic processes in the equation of motion 
of the sigma field, is given by the interaction with the heat bath. It depends on the temperature of the heat 
bath and consistently on the sigma field itself via the masses of the quarks and of the sigma mesons. 

• To account for the energy-momentum exchange between the two sectors we constructed an energy- 
momentum tensor, which is conserved for the full 2PI effective action. The approximations used for the derivation 
of the dynamics of the sigma field generate additional terms in the energy-momentum balance. 

In the present work we implement this consistent coupling numerically and study the full nonequilibrium dynamics 
of the entire system with respect to effects of the phase transition. A detailed description of the theoretical aspects 
of the model can be found in I29I1. 



1 The last term in the l.h.s. of equation ((TJ is written in the rest frame of the fluid. To restore a covariant form of the equation one should 
replace dt — ► u^d^a, where m m is the 4— velocity of the fluid element. 
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B. The linear sigma model with constituent quarks 



We use the linear sigma model with constituent quarks as a low energy effective model of QCD. It exhibits some 
important features of the suggested chiral phase structure of QCD. It has a crossover transition at vanishing (i and a 
first order transition line at high /i and lower temperatures, which terminates in a critical point at T c ~ 99 MeV and 
fx c ~ 207 MeV [39]. The quantitative parameters of the phase transition in the linear sigma model may differ from 
the values for the true QCD. Therefore, we expect only a qualitative similarity to the true QCD phase transition. 

The Lagrangian of the theory reads 



C =q [ij^d^ - g{a + 1757-7?)] q 



(3) 



Here q = (u, d) is the constituent quark field, a and 7? are the sigma and the pion fields. The strength of the coupling 
between the quarks and the chiral fields is g. In the vertex for the pion-quark coupling the 75 matrix enters to account 
for the pseudoscalar nature of the 7r-mesons and the isospin matrix f accounts for the iso- vector nature of the pions. 
The interaction between the chiral fields is given by the potential 



U(a,n) = ^ (a 2 



i^ 2 ) 2 — h q cr 



U . 



(4) 



Usually the parameters in (j4j are chosen such that the vacuum properties of the system are reproduced. Chiral 
symmetry is spontaneously broken in the vacuum and the sigma field acquires a finite expectation value of (a) = 
fix = 93 MeV while the pions have zero vacuum expectation value (7?) = 0. The explicit symmetry breaking term 
is h q (T — f^m^ with = 138 MeV. Thus v 2 = / 2 — to7t 2 /A 2 . Choosing A 2 = 20 yields a sigma mass m a = 



V2A 2 / 2 



ml 



604 MeV, which is within the mass range assumed for the sigma state in the meson spectrum. In 



order to have zero potential energy in the ground state Uq = m^/(4A 2 ) — / 2 m 2 - At a coupling g = 3.3 the constituent 
quark masses in vacuum are m q — 306.9 MeV. 

In the following we will use for the parameters the values given above - except for the coupling constant g. In 
principle, one can introduce a finite baryo-chemical potential and vary it to g et either a crossover, critical point or 
first-order transition. For our qualitative study, however, we follow (40l . |41| and vary the coupling strength g to 
consider the different possibilities. 

The grand-canonical partition function of the quarks in presence of the chiral field is 



Z = I VqVqVaVif exp 



d 4 xC 



(5) 



In the mean-field approximation Z can explicitly be calculated. The thermodynamic potential to one-loop level is 
given by 



V eS (a, 7?,T) 



T 
V 



InZ = U (a, 7?) + flq q (T, a, 7?) 



, (6) 



where the degeneracy factor d q = 12 is the product of Nf — 2 flavors, N c = 3 colors and two spin states. The energy 
of the quarks 



E q = Jp 2 



mi 



J q — y f ■ '-q 

depends on the chiral fields via the dynamically generated quark masses 

m = g (a + tt ) . 



(7) 



(8) 



At [i = the strength of the phase transition can be tuned by varying the strength of the coupling g. For g = 3.3 
the transition is a crossover, while for g = 3.63 the potential becomes flat at T c = 139.88 MeV as for a continuous 
transition (critical point). And for g — 5.5 the phase transition is discontinuous (first order phase transition). Here, 
one finds two degenerate minima at T c — 123.27 MeV corresponding to the two coexisting phases in a first order 
phase transition. 
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For the fluid dynamical description of a phase transition we need an equation of state, which relates the pressure to 
the energy density. Since the sigma field is not necessarily at its equilibrium value the equation of state will depend 
explicitly on the local value of the sigma field <r(x). This is different to standard chiral equations of state used in fluid 
dynamic calculations. The pressure of the quarks is locally given by 

p(a,7r,T) = -V e g(a,7r,T) + U(a,7r) . (9) 

The equation of state is completed by obtaining the local energy density from the thermodynamic relation 

e(a,ir,T)=T — p(a,ir,T) . (10) 

In the linear sigma model the sigma field is the order parameter of the chiral phase transition. In the rest of this 
paper we will neglect the fluctuations of the pion field and set it = (n) = for all times. 



C. Inclusion of dissipation and noise 

The damping coefficient r\ in the Langevin equation for the sigma mean-field was calculated in [29J for the linear 
sigma model with constituent quarks. Since we are mainly interested in the long-range fluctuations of the sigma 
modes we use the damping coefficient for the |fc| = mode in the Langevin equation. It is 



IT 



1 - 2n T 
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1 



>>>?, 



4 



3/2 



(11) 



It depends on the local value of the sigma field via the sigma mass and the quark mass. We calculate these masses 
by using the local equilibrium value of the sigma field a cq . It is obtained by minimizing the effective potential V s 



dV cS 



= and 



d 2 V cS 



d 2 c 



> 0. 



The sigma mass is given by the curvature of the effective potential at the minimum 



XT) 



d 2 V eS (T,a) 



da 2 



(12) 



(13) 



Alternatively to the use of the equilibrium value one could try to use the actual value of the sigma field for the 
evaluation of the quark mass in ([8]). For the determination of the sigma mass one presumably would have to study 
the (local) response of the sigma field to small variations on the basis of the equation of motion ((T|) . At the critical 
point the equilibrium sigma mass becomes very small, but we expect to see fluctuations, which can locally change this 
effective mass. Concerning a first-order phase transition it is known that in nonequilibrium in the spinodal region 
the sigma mass square becomes negative. In principle, the explicit cr-dependence of 77 = T](a, T) could take these 
nonequilibrium effects into account. Then, several subtle questions need to be addressed, e.g., which "mass" should 
enter formula (jllj) . if the effective mass becomes tachyonic. In the present work, however, we disregard these issues 
and use the local equilibrium value for the determination of the sigma and quark masses which in turn enter the 
damping rate (fTT|) . 

This results in a temperature dependence of the damping coefficient r\ as shown in figure [TJ At lower temperatures 
it vanishes because the decay of the zero mode of the sigma field in a quark-antiquark pair is kinematically forbidden 
as m a (T) < 2m q (T) — 2ga cq (T). Albeit confinement is not included in the model, the vanishing of the damping 
process a — > qq at low temperatures is i n ag reement with confinement. Physically the sigmas get further damped by 
interactions with the hard chiral modes [171 . |20| , which are not included in our mean-field approach. To account for 
them we add to 77 a value 2.2/fm [T2|. This contribution is attributed to the reactions a ■<-» 2tt. Consequently this 
damping effect should vanish when the sigma mass becomes smaller than 2m T . This occurs only in a critical point 
scenario, where the sigma mass is very small. We note also that?7 is very large, especially for the investigated scenario 
with a first order phase transition. This is due to the perturbative character of the derivation of rj and is manifested 
in the g 2 dependence of the result. This was discussed in [29| and is considered to be improved in future work. 

To summarize, for m a {T) > 2m q (T) — 2ga cq (T) the damping coefficient is given by equation (fTTj) . For 2m q > 
m a (T) > 2m^ we use the value 77 = 2.2/fm and for m a (T) < 2m w ,2m 9 r\ = 0. The stochastic field in the Langevin 
equation (flj has a vanishing expectation value 



(£(*)>« = o, 



(14) 



5 



first order phase transition 



critical point 

crossover 



Ill) 120 130 140 150 160 170 180 190 200 
T [MeV] 



FIG. 1: Temperature dependence of the damping coefficient 77 for a the different couplings g = 5.5, g — 3.63 and g — 3.3, which 
correspond to scenarios with a first order phase transition, a critical point and a crossover [29(. 



and the noise correlation is given by the dissipation-fluctuation theorem 



(15) 



D. Source term for the quark fluid 

In [29j we constructed a conserved energy-momentum tensor of the entire system using the formalism of the 2PI 
effective action. For any approximation to the full quark propagator the situation is more difficult due to the space- 
time dependence of the effective mass generated by the dynamic symmetry breaking. The total energy-momentum 
tensor can be represented as a sum of two contributions 

where the first term comes from the quark-antiquark fluid and the second term from the sigma field. 
The energy-momentum conservation equations, d^T^ = can be written as 

d^{x) = - WOr) = S»(x) , (17) 

where S" is the source term for the fluid dynamics. In [2fj| we have obtained the following approximate expression 
for the part of the sigma field 

d^Zvpri*) = i-9Ps(x) - r){x)d t a{x)) d v a{x) . (18) 

It consists of the standard mean-field result including the scalar density p s and the correction term given by the 
damping coefficient. It can, however, not account for the average energy transfer from the heat bath to the field given 
by the auxiliary noise field £. This issue was discussed in [42j]. 

The quark contribution can be presented in symmetric form in terms of the quark propagator (dependent on center 
and relative variables x and u of the arguments) 

wSppt.Os) = -fa(d:s: h (x,u)\ u = Q r+dzs: h (x,u)\ u=01 » + ... (19) 

The first term can be evaluated explicitly, 



T q W = ~\ (d" u S: h -(x,u)\ u = r + dZS: h (x,u)\ u = 0l n 



2 



(20) 



d 3 p p^p" 
(2tt) 3 p" 

It gives the energy-momentum tensor for an ideal fluid with the energy density and the pressure obtained from the 
equilibrium one-loop effective potential in mean-field approximation. This is exactly what we intend to use for the 
fluid dynamic description of the quark-antiquark fluid. Higher-order corrections to the energy-momentum tensor of 
the quark fluid are neglected in the present set up. Then, the source term in equation ^ can be identified with 

= — ^/i-^appr. • (21) 
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III. NUMERICAL IMPLEMENTATION 



A. Fluid dynamics 

In the relativistic fluid dynamic equations ([2]) the energy-momentum tensor of an ideal fluid reads 

T" v = (e + p)u»u v - pg^ , (22) 

where e is the local rest frame energy density, p is the local rest frame pressure, = 7(1, v) is the local four- velocity 
of the fluid, 7 = (1 — v 2 )^ 1 / 2 and g^" — diag(l, —1, —1, —1) the metric tensor. 

For the solution of © we use the full (3+l)d SHarp And Smooth Transport Algorithm (SHASTA) ideal fluid 
dynamic code [43], Ejj . The causal transport of matter is assured on a numerical level by fulfilling the Courant- 
Friedrichs-Levy criterion At/ Ax = A < 1 and the SHASTA code requires values A < 1/2 for numerical stability. We 
use A = 0.4, Ax = 0.2 fm and thus At = 0.08. Our grid volume is (128Ax) 3 . 

The numerical implementation of the source term (|21[) is as follows. After performing the fluid dynamic step 
in the standard fashion for S u — we subtract the source term S v from the energy and momentum density in the 
computational frame. Then, we can again calculate the local rest frame quantities. This gives a very good conservation 
of energy and momentum of the entire system as we will show in section fill CI The latter method also proved to work 
well in multi- fluid dynamics (45j . 

For the equation of state we do not obtain a simple relation p(e) in the nonequilibrium evolution because the sigma 
field is not fixed at its equilibrium value. Therefore, the pressure ^ and the energy density (|10p depend explicitly on 
the local value of the field, which can be viewed as an external parameter in the thermodynamic sense. Technically 
the following has to be done: With the energy density from the fluid dynamic calculation at a given point a;, efl u ;d(a;) 
the equation for the thermodynamic energy density (| L0[) needs to be inverted taking into account the local value of 
the sigma field a(x). The local temperature T(x) is given by the solution of 

efluid(z) ~ e(cr,T) = , (23) 

and used to calculate the thermodynamic pressure ([9]) . For the transformations between the local rest frame of a fluid 
cell and the computational frame the equation of state is accessed very often in each time step. Since it is very time 
consuming to invert (fTU|) numerically, we use a properly parametrized equation of state. 



B. Initial conditions 



For a qualitative investigation the initial conditions are kept simple. We choose an initial temperature profile, 
which is uniform in z-direction over a length l z = 6 fm and ellipsoidal in the x/y-plane. The ellipsoidal shape should 



mimic the overlap region in a heavy-ion collision. Its major and minor axes are b = \Jr\ — b 2 /A and a — ta~ b/2, 

where b = 6 fm is the supposed impact parameter and ta = 6.5 fm the radius of Au nuclei. The temperature is 
smoothly distributed over this ellipsoidal region by a Wood-Saxon like distribution with the maximum at the initial 
temperature T; n ; = 160 MeV, which is well above either of the phase transitions considered here, 



T{x,t = 0) = 



(1 + exp((r - R)/a)){l + exp ((\z\ - l z )/a)) 



(24) 



with a surface thickness of a = 0.3 fm, f = \/ x 2 + y 2 and 



abr 



, for f 7^ 

f> = ) v /b 2 x 2 +a 2 y 2 r 

a for r = 



(25) 



By minimizing the effective potential the equilibrium value of the sigma field a cq is found. The thermal equilibrium 
state has Gaussian fluctuations around the expectation value, the variance of which is 



^ 2 ) = l\- 
V m% 

Thus, the sigma field is initially Gaussian distributed around its equilibrium value, 



P(a(x)t = 0) = 



1 



y/2ir(5<r 2 



: exp 



(a(x) - Cr eq (a;)) 2 



2(Sa 2 



(26) 



(27) 
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According to (flOl) the initial energy density of the quark fluid can be calculated with ((24)) and (|27|) . 

The velocity profile is v z (x,t = 0) = \z\/l x ■ v max , with w max = 0.2, while initially there are no transverse velocities 
v x (x,t = 0) = v y (x,t — 0) = 0. Before we start the fluid dynamic expansion the equation of motion for the sigma 
field is solved a couple of times at Ti n j to generate an initial distribution of the time derivative of the sigma field. The 
pion field is initially set to zero and neglected during the simulation, since we expect that mainly the sigma field as 
the order parameter of chiral symmetry is affected by the phase transition. 



C. Energy- momentum conservation 

The energy-momentum conservation was discussed in [29j . where we found a conserved energy-momentum tensor 
for the full 2PI effective action. In the truncated theory additional terms contribute to the energy-momentum balance. 
We first investigate the overall energy-momentum conservation before turning to the evolution of the system. 

The correction to the divergence of the energy-momentum tensor of the sigma field in (|18p is given by the dissipative 
part of the equation of motion. It describes the energy transfer from the field to the fluid caused by the interaction 
between the fields and the fluid given by the pseudoscalar density p s and by dissipation for each numerical time step 
At 

AE n = (gps(x) + T](x)d t a(x)) d t a(x)At . (28) 

In addition, the presence of the stochastic field £ causes an average energy transfer AE% from the heat bath to the 
field. We determine AE^ by comparing AE^ to the numerically obtained energy difference in the sigma field before 
and after each time step of the simulation. 
The total energy of the sigma field is given by 

E( , = l(d t a) 2 + ±(Vaf + U(a). (29) 

It has a kinetic, spatial fluctuation and potential energy term. Then, the following expression is used as an energy 
source term in the numerical simulation 

S° +i = ^ (AE V + AEs) . (30) 

In [42j we studied the relaxation of the sigma field in a static heat bath including this source term to account for 
the energy exchange. Here we present calculations for the expanding fluid where the momentum conservation must 
be taken into account, too. The spatial part of the source term, which is responsible for the momentum exchange, is 
calculated in a similar way. The momentum transfer due to the interaction and dissipation is given by 

AM, - (gp s (x)+f 1 (x)d t a(x))d<j(x)At, (31) 

and AM^ is obtained from the comparision to the change in the momentum density of the field 

M a = d t odo . (32) 

Then, the spatial part of the source term reads 

Sr, +£ = ^ (AM, + AM { ) . (33) 

On the fluid dynamic side of the coupled evolution we apply the approximation of an ideal fluid by neglecting any 
corrections to the energy- momentum tensor of the quarks given by equation (|22[) . 

In the beginning of the simulation the field configuration gains energy due to increasing fluctuations. In some cells 
this increase in field energy may exceed the energy in the fluid. In these cells the energy density is set to zero at the 
end of each numerical time step. In the beginning we, thus, observe a net increase of the total energy by < 10% in the 
scenario with a first order phase transition, while the energy conservation is well fulfilled during the expansion. In a 
scenario with a critical point the energy is very well conserved during the entire evolution. The quark fluid reaches the 
edges of the grid at around t ~ 8 fm and disappears. In a test case with a larger grid the total energy was conserved 
well for a longer time. In addition, the evolution of the system in the inner region was not altered. The results are 
shown in figure [U 

In addition to the energy conservation we also checked the momentum conservation. The total momentum of the 
system is close to zero. Due to the finite initial time derivative of the sigma field each initial field configuration has a 
small overall momentum, which is slightly enhanced in the beginning of the simulation. The total momentum in each 
of the three directions is of the order of 0.1% of the total momentum in positive direction. 
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FIG. 2: The total energy of the system for a scenario with a first order phase transition (a) and with a critical point (b) It is 
the sum of the fluid energy in the laboratory frame and the field energy given in equation (|29p . 




FIG. 3: Time evolution of the fluid energy density along the a;-direction in the plane y = z 
point (a) and with a first order phase transition in (b) 



for a scenario with a critical 



IV. NONEQUILIBRIUM EFFECTS IN CHIRAL PHASE TRANSITIONS 

With a constant damping coefficient r\ one can generally investigates dissipation and relaxation 46] . In this section 
we investigate the evolution of the entire system under the temperature-dependent damping coefficient r\ — i](T). The 
temperature dependence of rj — rj(T) includes more aspects of the phase transition scenario than a constant damping 
coefficient. 



A. Time evolution of the coupled system 

The time evolution of the energy density and the sigma field in ^-direction and at y = z = is shown in figures [3] 
and [4] for both phase transition scenarios. In figures [5] and [6] we additionally show the longitudinal expansion along 
the z-direction and at x — y — 0. We see how the quark fluid expands and the system dilutes. For a first order phase 
transition we can observe the energy transfer from the field to the fluid in figure [3(b)] At x ~ and t ~ 9fm the sigma 



field finally relaxes, see figure |4(b)| At the same location and place we observe an increase in the fluid energy density 
in figure [3(b)] This effect is less pronounced in a critical point scenario. Here, relaxation of the sigma field occurs 
a lot faster, see figure 3(a)| and shows a qualitatively different behavior than for a scenario with a first order phase 



transition. During the relaxational process there are oscillations for the critical point scenario, because the damping 

coefficient vanishes for a certain region around the phase transition. 

The time evolution of the average value of the sigma field (a) and its fluctuations J (5a 2 ) = -J ((a — (er)) 2 ) are 
shown in figure [7] Figure E] shows the time evolution of the average temperature (T) and its fluctuations (ST 2 ) = 
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FIG. 4: Time evolution of the sigma field along the :r-direction in the plane y = z = for a scenario with a critical point (a 



and with a first order phase transition in (b) 




FIG. 5: T ime evolution of the fluid energy density alo ng the z-direction in the plane x = y — for a scenario with a critical 
point (a) and with a first order phase transition in (b) 



\J ((T — (T)) 2 ) . The average is taken over an initially hot and dense sphere with radius r = 3 fm. We first discuss 
the critical point scenario and then the scenario with a first order phase transition. 

The phase transition temperature in a critical point scenario T c = 139.88 MeV is crossed at around t = 4 fm. 
Fluctuations of the sigma field are largest between t = 2 fm and t = 4 fm and stay enhanced until around t = 10 fm. 
The average sigma field smoothly relaxes towards its vacuum value. Strong oscillations of the average sigma field occur 
during the relaxational process, when large parts of the system do not experience damping. The effective potential 
with a critical point is very flat at the transition temperature and reheating is not observed. Instead the cooling is 
slightly decelerated between t = 5 fm and t = 6 fm as seen in figure [SJ 

There are two reasons why the average sigma field in the scenario with a first order phase transition relaxes at later 
times than for a critical point scenario. First, the phase transition temperature of the first order phase transition, 
T c = 123.27 MeV, is lower than at a critical point. And second, the sigma field is more strongly damped at high 
temperatures in a scenario with a first order phase transition than in the critical point scenario. At the phase transition 
temperature T c the sigma mass drops below the threshold of quark-antiquark production but is larger than twice the 
pion mass. The damping coefficient shows a discontinuity from r/(T > Tc) ~ 26.6/fm to rj(T < Tc) — 2.2/fm. In this 
scenario the average of the sigma field stays constant up to almost t ~ 5 f m and the relaxational process starts only 
when the average temperature is below the phase transition temperature as can be seen by comparing figure [7] with [5] 
The cooling of the system is inhomogeneous. While the center is still very hot, outer parts have cooled down already. 
The fluctuations of temperature within the inner region with radius r = 3 fm have a variance of almost 10 MeV. 
We conclude that even in this small inner region the outer parts have already cooled below the phase transition 
temperatures at t = 5 fm and that these parts relax first due to the lower damping coefficient r\. When large parts of 
the system experience the smaller damping of r\ = 2.2/fm the average sigma field shows oscillations during the final 
relaxational process also for a scenario with a first order phase transition. 
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FIG. 6: Time evolution of the sigma field along the z-direction in the plane x = y = for a scenario with a critical point (a) 



and with a first order phase transition in (b) 



/ critical point, \J (Sa 2 ) 
first order, \J (Sa 2 ) 
critical point, (a) 
first order, (a) 




FIG. 7: The average values and the variances of the fluctuations y/ {5a 2 } of the sigma field for r/ — r/(T) and both phase 
transition scenarios. 



B. Supercooling and Reheating 



When the first order phase transition temperature is reached after t ~ 5 fm large parts of the system are still in 
the chirally symmetric phase as the average value of the sigma field is still (a) < 10 MeV. These large deviations of 
the sigma field from its equilibrium value is the nonequilibrium effect of supercooling. Due to the barrier separating 
the two minima of the thermodynamic potential, the sigma field is traped in the chirally symmetric state even at 
T < T c , until the barrier disappears. Then the sigma field rolls down into a lower minimum corresponding to the 





critical point, \J (ST 2 ) 






first order, y/ (<5T 2 ) 
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FIG. 8: The average values and the variances of the fluctuations a/ (ST 2 ) of the temperature for r\ = T](T) and both phase 
transition scenarios. 
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FIG. 9: The time evolution of the intensity of sigma fluctuations for a critical point scenario (a) and for a scenario with a first 
order phase transition (b) 



broken phase. 

The released potential energy is transformed effectively into kinetic energy of field oscillations, which leads to the 
dissipation of energy into the fluid via r](d t <j) 2 in the source term (fig]) . In figure[5]we can clearly observe the reheating 
effect at the first order phase transition. Between t — 7 fm and t = 9 fm the system in the central region is reheated 
from T ~ 114 MeV below T c to T ~ 124 MeV slightly above T c , followed by a subsequent cooling. Thus, the reheating 
causes the system to cross the phase transition boundary two more times, once in the reverse direction from the low 
temperature phase to the high temperature phase around t — 8 fm and again from above to below at around t = 9 fm. 
This explains the delayed relaxation of the average sigma field in the scenario with a first order phase transition 
compared to the critical point scenario. Moreover, during the relaxation and reheating process the fluctuations in the 
sigma field are enhanced between t — 5 fm and t = 8 fm. 



C. The intensity of sigma fluctuations 



In various studies on the formation of DCC a large amplification of pion fields was observed at the phase transition 
|22J. The pionic excitations were found to be triggered by violent oscillations of the sigma modes. Moreover, 
the zero mode of the sigma field is the order parameter of chiral symmetry breaking. The soft modes are thus expected 
to show large fluctuations at a critical point in thermodynamic systems. We are, therefore, especially interested in 
the intensity of sigma fluctuations at the phase transition. 

The intensity of the sigma fluctuations is given by (2|| |47| 



a k a k 



d 3 k (2%) 3 2u k (2ir) 3 2ujk 



(ul\5a k \ 2 + \d t o k f 



(34) 



where a* k and a k are the Fourier coefficients in the expansion of the sigma fluctuations around the equilibrium value, 
Sa(x) = <t(x) — cr eq (a;), and its conjugate momentum field dt<r(x). 

At later times when the nonlinearities in the equation of motion can be neglected the quantity N a in (|34[) gives the 
number of sigma particles produced from the excitations of the sigma field. The energy co k of the fcth mode of the 
sigma field is given by 



^fe = V\ k \ 2 + ™% 



(35) 



We define the sigma mass via the curvature of the effective potential at its equilibrium value corresponding to the 
average temperature T av in a hot region of radius r = 3 fm. It turns out that the intensity of sigma fluctuations does 
not depend strongly on the choice of the radius, as long as it does not extend too far into the surrounding vacuum. 

At the end of each time step the Fourier transform of the sigma field is calculated by a Fast Fourier Transform 
algorithm. These algorithms are implemented effectively only on lattices with N = 2™, n € N, cells in each direction. 
This limits the realistic choice of how to define the entire volume of the system. We perform Fourier transformations 
on the entire grid. 

For the two phase transition scenarios we show the intensity of sigma fluctuations for the zero mode and low 
momentum modes in momentum bins of A|fc| = 50 MeV. Finally, we compare these values to the development of 
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t [fin] 

FIG. 10: The time evolution of the deviation of the sigma field from its equilibrium value for a critical point scenario and a 
scenario with a first order phase transition. 



the deviations of the sigma field from its thermal equilibrium value averaged over the initially hot and dense region. 
The initial conditions, see section 1111 B( are chosen such that the sigma field is in equilibrium with the quark fluid 
at the initial temperature distributed according to (|24|) . In both phase transition scenarios the intensity of sigma 
fluctuations is increased during the expansion and cooling see figure [HI For a critical point scenario the results are 



shown in figure 9(a) The intensity shows a peak at t = 2 fm and a second but smaller peak at t = 6 fm. Between 
t = 2 fm and t = 6 fm the deviations of the sigma field from its equilibrium value are enhanced, see figure 1101 

For the scenario with a first order phase transition the time evolution of the intensity of sigma fluctuations is shown 
in figure [9(b)] Figure [TU] shows the deviations of the sigma field from its equilibrium value. Since these deviations are 
larger than for the scenario with a critical point, the intensity of sigma fluctuations is also larger. Largest deviations 
are seen around t — 6 fm, where large parts of the system are supercooled. Then, the system quickly relaxes to a new 
ground state. This process is accompanied by reheating of the heat bath. This explains a shoulder, before the system 
relaxes and the deviations vanish. 



V. SUMMARY AND OUTLOOK 



We presented a consistent dynamic model to describe nonequilibrium phase transitions in heavy-ion collisions. 
We included damping and noise terms originating from the interaction of the chiral fields with the quark fluid. To 
ensure energy and momentum conservation of the coupled system we introduced a source term into the equations 
of relativistic fluid dynamics and thus explicitly took the back reaction on the heat bath into account. The various 
aspects of energy and momentum conservation were discussed in detail. We investigated the expansion of the coupled 
system with a damping coefficient rj — rj{T) for a scenario with a critical point and with a first order phase transition. 
Nonequilibrium effects like supercooling and reheating are clearly observed in case of the first order phase transition. 
The number of coherently produced soft sigmas is much larger in a scenario with a first order phase transition than 
in a scenario with a critical point. 

Several extensions of the presented framework are conceivable which would make the approach quantitatively more 
realistic. Most straightforward is the inclusion of soft pions as classical fields subject to stochastic dynamics and of 
hard pions which contribute to the heat bath. Possible extensions include gluonic degrees of freedom in the form 
of a dilaton field [22[ or the Polyakov loop [48]. This leads presumably to a more realistic phase diagram on the 
temperature-chemical potential plane. Finally a finite baryo-chemical potential can be introduced which allows to 
access different types of the phase transition without varying the coupling constant. The comparison of the fluctuation 
dynamics in the static and dynamical background is a particularly interesting topic for future studies |27j | . 
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